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Abstract. The behaviour of extended particles with exclusion interaction on a 
one-dimensional lattice is investigated. The basic model is called £-ASEP as a 
generalization of the asymmetric exclusion process (ASEP) to particles of arbitrary 
length £. Stationary and dynamical properties of the £-ASEP with periodic boundary 
conditions are derived in the hydrodynamic limit from microscopic properties of the 
underlying stochastic many-body system. In particular, the hydrodynamic equation 
for the local density evolution and the time-dependent diffusion constant of a tracer 
particle are calculated. As a fundamental algebraic property of the symmetric exclusion 
process (SEP) the S , J7(2)-symmetry is generalized to the case of extended particles. 
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1. Introduction 

There is renewed interest in the investigation of extended particles with exclusion 
interaction. The basic model, which will be referred to as the £-ASEP in the following, is 
a generalization of the well-studied asymmetric simple exclusion process (ASEP) [1,2]. 
It describes the motion of hard rods in one-dimensional discrete space by extended 
particles which move along a lattice according to stochastic hopping dynamics. 

The original concept of the £-ASEP was introduced in 1968 in a paper by 
MacDonald and Gibbs treating protein synthesis [9,10]. During this process, ribosomes 
move from codon to codon along a m-RNA template, reading off genetic information and 
thereby generating the protein step by step. The ribosomes are modeled as extended 
particles, which hop stochastically along a chain without overlapping each other. Each 
particle covers several adjacent lattice sites to account for the blocking of several codons 
by a single ribosome. The attachment of the ribosomes to the m-RNA for the initiation 
of the protein synthesis and their detachment at the point of termination are modeled 
by open boundaries, where particles may enter and exit the lattice with rates that differ 
from the bulk hopping rates. Using mean-field theory, the authors studied the steady 
state of this process. More recently the time-dependent conditional probabilities [3], 
the dynamical exponent [4] and the phase diagram of the open system have been 
determined [5-8]. 

However, the understanding of symmetries of the model and of its hydrodynamic 
limit has remained incomplete so far. In [6] a hydrodynamic equation is proposed 
phenomenologically, employing fitting parameters, which are matched to simulation 
data. In the present paper, several basic physical and mathematical properties 
of the £-ASEP with periodic boundary conditions are derived from its microscopic 
dynamics, generalizing a mapping to the zero range process [11] and employing quantum 
Hamiltonian techniques [2]. In particular, we obtain the hydrodynamic limit governing 
the density evolution of the £-ASEP on the Euler scale. 

The outline of this work is as follows: After introducing the two fundamental 
models (£-ASEP and zero range process) and reviewing some facts about their stationary 
properties in section 2, the dynamics of the £-ASEP are studied by two different 
approaches: In section 3, the investigation of the motion of a tagged particle in the 
framework of the quantum Hamiltonian formalism leads to an expression for the average 
velocity and the time-dependent diffusion constant of the tracer particle. The velocity 
term is then confirmed by the general form of a hydrodynamic equation of the £-ASEP 
which is derived in section 4. Section 5 finally exposes the hidden 5'[/(2)-symmetry as a 
fundamental algebraic property of the £-SEP, which arises from the £-ASEP by requiring 
left /right-symmetric hopping rates. 
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Figure 1. Two approaches to model the £-ASEP (left: extended particles which cover 
I lattice sites each and move via next-neighbour hopping, right: point particles with 
certain initial states and modified hopping rules.) 



2. £-ASEP and ZRP 

2.1. Thei-ASEP 

The £-ASEP is the discrete nonequilibrium analogue of a one-dimensional Tonks gas, 
including the ASEP as a special case for £ = 1. iV particles are placed on a one- 
dimensional lattice S consisting of L sites k = 1 . . . L (Fig. 1(a)). Each particle covers I 
adjacent sites. The parameter i is an integer number which determines the extension of 
the particles in units of the lattice spacing a. In the following, t will be called the length 
of a particle. As time proceeds, the particles change their locations on the lattice by 
next-neighbour random hopping under exclusion interaction. Provided that their right 
and left neighbour sites respectively are not occupied, they move one site to the right 
with rate p, or one site to the left with rate q. The location of a particle on the chain 
is denoted by the location of its right end. 

In figure 1 this model is compared to an equivalent description (lb), where each 
extended particle is composed of i monomers. In the latter case, the initial states are 
restricted to those where particles are grouped in ^-tuples at adjacent sites. Hopping 
processes take place between sites which are i lattice spacings afar and which enclose 
I — 1 occupied sites in between. This monomer description is useful for the application 
of the quantum Hamiltonian formalism to the £-ASEP (see below). 

2.2. Mapping between the I- ASEP and a ZRP 

One important steady-state property of the ^-ASEP follows directly from its definition. 
As sites are always correlated, the £-ASEP does not possess a stationary product 
measure. However, the existence of a stationary product measure is an important 
ingredient in the conventional derivation of hydrodynamic properties. In order to recover 
it, the £-ASEP can be mapped onto a different lattice gas model: the zero range process 
(ZRP) [12]. The ZRP does have a stationary product measure, a fact which is used in 
section 4 to derive the hydrodynamic equation for the density evolution of the £-ASEP. 
Furthermore the ZRP-picture will be of help when considering the motion of a tagged 
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Figure 2. Mapping between ZRP and £-ASEP: Lattice sites are turned into particles, 
particles are replaced by holes; the stochastic hopping rates p and q are interchanged. 




particle in section 3. The zero range process is named after the fact that its particles 
have zero interaction range, i.e., there is no exclusion and jump rates do not depend on 
the occupation number of the target site. 

In the following, a ZRP will be considered where particle hopping from a site 
occupied by n particles occurs with fixed biased rates q and p to the left or right 
respectively. Like in the case £ = 1 [11], the £-ASEP can be mapped onto this ZRP by 
replacing particles by ZRP sites and holes by ZRP particles (Fig. 2). The appropriate 
coordinate transformation (compare figure 3) between the ZRP with M particles on a 
lattice of N + 1 sites j — 0, 1, . . . , N, and the £-ASEP, having N particles on a lattice 
of L = £N + M sites k — 1, . . . , L is explicitly given by: 

+ j£ (la) 

where k is the £-ASEP lattice site corresponding to a certain ZRP site j and Cj(t) 
denotes the ZRP particle density at site j at time t. The ZRP site j = is not turned 
into a particle but determines the position of the first £-ASEP particle. This procedure 
guarantees the uniqueness of the transformation. In the continuum limit, where the 
lattice constant a approaches zero, the discrete coordinates j and k may be replaced by 
continuous variables y and x: 

x= f V duc(u,t) + ly + ^[c(0,t)-c(y,t)]+O(a 2 ). (16) 
Jo I 

The ZRP-density c is related to the £-ASEP particle density p by: 
and 

c(y) + £ 

Whenever speaking of the density p in the following, indeed the particle density, i.e. the 
fraction of particles per lattice unit, is referred to, as opposed to the coverage density 
p c and hole density p h respectively. 
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Figure 3. Transformation of coordinates and densities: ZRP and £-ASEP quantities 
are related by a transformation between time-independent and time-dependent 
coordinates; the ZRP-site 'zero' provides a left boundary of the £-ASEP lattice (shown 
here for N = 5, M = 6, 1 = 2). 



2.3. Stationary state 

The stationary state of a zero range process is known [12]. For a periodic system all 
ZRP configurations of the present model are equally probable. Due to the existence of 
the one-to-one mapping between ZRP and £-ASEP configurations on (periodic) lattices 
of fixed length and particle number (N, L — IN) and (L, N) respectively, the stationary 
weights of the £-ASEP must also be distributed equally among all configurations on a 
ring. 

The stationary properties of such a system can be deduced from a partition function 
of the form 

ax 

Z=Y, * N Zn (3) 

where z denotes the fugacity and is the A^-particle partition sum. N max indicates 
the maximum number of particles fitting completely on L sites. As all states contribute 
equally, Zn is given by the number of possible different A^-particle configurations on a 
lattice of length L [5, 13]: 

H l - ( ;- i)jv ). w 

The expectation value and the fluctuations of the particle number are calculated from 
the first and second derivatives of Z with respect to the fugacity z using the following 
standard relations: 

(AO = Lp = N / Zn = z^lnZ (5a) 

Z dz 

(N*) - (AT) 2 = 4 (z±lnz) . (56) 
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In the hydrodynamic limit the partition sum is approximated by its maximum term and 
a stationary density-fugacity relation may be derived as: 
p(l-(l-l)p)<-i 

z = — (i^Ipy — • (6) 

The stationary density fluctuations are given by: 

For £ — 1, equation (7) reduces to the well-known expression for the compressibility of 
the ASEP 

kasep = (iVW ~ {N)2aSEP = P(l - P) ■ (8) 

The extra factor in (7) may be written as j(L — {£ — 1)N). This term, ranging between 
l/£ and 1, specifies which fraction of the system is formed by holes and particle ends. It 
accounts for the fact that the extended particles, constructed as a chain of £ monomers, 
are stiff, and that they move simultaneously. The fluctuations in the particle number 
taken per volume of holes and ends (L — (£ — 1)N) are the same for any monodisperse 
system which contains particles of an arbitrary length £^ provided that the particle 
density and the hole density are fixed a priori. 



3. Motion of a tagged £-ASEP particle 

In experiments, a common procedure to investigate the dynamics of a diffusive system is 
to mark a particle and to track its motion. In the following, the motion of such a tagged 
particle is examined with analytical tools. The simplest conceivable case is the one of a 
system without interparticle interactions and without the influence of an external field. 
A free particle moving according to symmetric hopping rates on a lattice performs a 
random walk. Starting at time t = at an initial position x , its location at a time 
t > fluctuates stochastically around the expectation value (x) = xq. An external field, 
i.e., biased hopping rates, causes a particle to move forward into the direction of the 
drive. The expectation (x — xq) of the distance covered by this particle is different from 
zero and proportional to its average velocity. Diffusive fluctuations arise again due to 
the stochastic motion. 

If one traces a certain particle of the £-ASEP, a similar type of motion is expected 
which may be decomposed into a drift term and a diffusion term. However, the 
calculation of the corresponding average velocity v and the diffusion constant D is not 
trivial anymore because collisions cause time-correlations in the random displacement of 
the tagged particle. In order to obtain the average velocity of the marked particle, one 
counts the hopping events that this particle accomplishes starting from a time t — 0. 
At a time t > 0, the distance (in units of the lattice spacing) the particle has covered 
is given by the difference of the number of left and right jumps it has performed. The 
average distance is proportional to the average velocity v of the tagged particle over 
that time. 
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In the ZRP-picture, the average velocity of £-ASEP particle number % corresponds 
to the difference of leftward and rightward ZRP currents j(i — > i — 1) — j(i — 1 — > i). 
The calculation of these currents is carried out conveniently in the quantum Hamiltonian 
Formalism where the master equation of the ZRP is rewritten in terms of a Schrodinger 
equation in imaginary time [2]: 

±\P(t)) = -H ZRP \P(t)) . (9) 

The ZRP Hamiltonian H ZRP acts on the state space H of the probability vectors \P(t)), 
which are linear combinations of the probabilities P v (t) to find the ZRP lattice in a state 
i] at a time t 

\p(t)) = Y,m\v)- (io) 

v 

The ZRP Hamiltonian takes a simple form in terms of single-site operators. Products 
of those operators represent transitions between different system configurations, i.e. 
hopping events between adjacent sites. Let aj and alf be operators which create and 
annihilate a particle at site j. The off-diagonal part of the Hamiltonian is formed by the 
negative sum of all such hopping operators, multiplied with the corresponding hopping 
rates q and p: 

N-l 

H ZRP,off = £ h off 

u°ff _ + - - + 

ij+i ~~ Q a j a j+i ~ P a j a j+i ■ 

The diagonal part is deduced from the condition that probability must be conserved, 

i.e. {s\H ZRP — , (s\ being the (N + l)-dimensional constant summation row vector 

(1, 1, . . . , 1). The relations 

(s\aj = (s\l t 

(12) 

\s\aj = {s\rrij 

where rrij is a diagonal operator replacing in its action on (s\ yield: 

N-l 

ffW= §^ +1 (13) 
h jtj+ i = -qajaj +1 - pajaj +1 + qrrij + pm j+1 . 

3.1. Extension of the state space 

Focusing again on the task to keep track of local hopping events, this Hamiltonian is 
modified as follows. In order to create a counting mechanism for particles hopping 
between the chosen sites i,i — 1, one enlarges the tensor state space 7i of the ZRP 
lattice by an additional infinite dimensional subspace K, whose basis vectors represent 
the number k of backward minus forward hops between % and % — 1. Their entries are 
all zero except the k th counted from the "middle" downwards: 

K = Span{|A;) : k G Z} (14) 
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One modifies the Hamiltonian such that its action on the elements of K transfers \k) in 
|/c=F 1) whenever its action on 7i leads to a hopping event between the lattice sites i and 
% — 1 in the ZRP. This is achieved by extra terms 



H = H 



ZRP 



+ qa~ i 



{l-x-)+pa i _ 1 aJ(l-x+). 



(16) 



(17) 



All ZRP operators are enlarged to the new state space such that their action remains 
local on 7i. x_ and x + are ladder operators which act nontrivially only on K, i.e., 
the zero range operators and the counting operators commute. x± are defined by the 
relations: 

x_\k) = \k-l) 
x+\k) = \k + 1) . 

Finally one constructs an operator x on K for which all \k) are eigenstates with 
eigenvalues k: 

x\k) = k\k). (18) 
In the chosen basis (15) of N, x is represented as a diagonal matrix: 

x = diag(...,-2,-l,0,l,2,...). (19) 

The average velocity v and the diffusion coefficient D are derived from the long time 
limits of the expectation value of x and its fluctuations: 



v = lim ^r(x) 

D = hm — 

t-*oo at 



The n th moment of the operator x in the stationary state \P*) is given by: 
(x n ) = (s\x n exp(-Ht)\P*,0) 



(20) 
(21) 

(22) 



where 



|P*,0) 



\P*) ® |0) 
(s\ = ( Sl \ ® (s 2 | 

( s ,h(4M-)e^, 

While the bra-vector in (22) is the summation vector, adapted to the new Hilbert 
space, the ket-vector is given by the tensor product of \P*) G 7Y\K, characterizing the 
stationary state of the ZRP, and |0) £ N, which sets the counting mechanism to zero at 
the beginning. 
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3.2. General calculation of moments 

In order to calculate the moments (x n ) of the distribution of the tagged particle 
explicitly, the new Hamiltonian (16) is formally split into a non-pertubative term H 
and a pertubative part V, where 

H = H ZRP (23) 

V = qa+_ ± a~(l - x_) + pa'^a^l - x+) . (24) 

The exponential exp (—Ht) is expanded in a time-ordered Dyson series: 

exp (-(H + V)t) = exp (-H t)[l - f dnVfr) (25) 

Jo 

+ /*dn r dr 2 \/(n)\/(r 2 ) - . . .] 
jo Jo 

where 

V(r) = exp (H r)Vexp (-H r) . (26) 

For the derivation of the n th moment of x, only terms which are at most of order n have 
to be taken into account. All higher orders vanish identically. Thus 

(x) = (s\xex.p(-Ht)\P*,0) 

= (s\xexp(-H t)\P*,0) 

-(s\xexp(-H t) I' drV(r)\P*,0) 
Jo 

= (x) (0) + (27a) 

and 

(x 2 ) = (s\x 2 exp(-Ht)\P*,0} 

= (s\x 2 exp(-H t)\P*,0) 

-(s\x 2 exp(-H t) C drV(r)\P*,0) 
Jo 

+ (s\x 2 exp(-H t) fdnVfa) r dT 2 V(r 2 )\P*,0} 
Jo Jo 

= (a?)<® + + (a?)® . (27b) 

Setting (x(t = 0)) = 0, these expressions reduce to 

(x) = (28a) 
(x 2 ) = (x 2 )^ + (x 2 )^ . (28b) 

Explicit expressions for the higher moments in terms of correlation functions can be 
obtained analogously, but are not considered below. As this strategy of calculating the 
moments of a random walk driven by some underlying Markov process does not appear 
to be widely known, we stress that this derivation can be applied to any counting process 
for any underlying Markov chain given by some H . The perturbation V is determined 
by the events that increase (multiplied by x + — 1) or decrease (multiplied by x_ — 1) the 
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counter (particle position) respectively. In particular, the generating function (exp (ax)) 
can be obtained by replacing x± — * exp (±a) in V. This proves that indeed the m th 
moment is given exactly by the terms up to same order in the perturbative expansion. 
Notice that up to here we did not use stationarity of the distribution P*. The expressions 
derived above are valid for any initial distribution. 



3. 3. Driven motion: Calculation of the average particle velocity v 

First an explicit expression for the average particle velocity (20) shall be obtained. We 
present the calculation in some detail. Equation (28 a) implies: 

(x) = -(six exp (-.Hot) [* <1tV(t)\P*,0) 

Jo 

= -(s\x f drV(r)\P*,0) (29) 



where exp (—Hot) has been absorbed into (s\. Substituting the time-dependent operator 
V(t) by (26) yields: 

(x) = -(six f dr exp (H t)V exp (-H t)\P* , 0) 
Jo 

= -(six I drV\P*,0) . (30) 
Jo 

The integration is then trivial: 

(six f drV\P*, 0) = t(s\xV\P*, 0) . (31) 
Jo 

Inserting the explicit form of V and applying relations (12) for the action of the creation 
and annihilation operators results in the simplified expression for the expectation value 

(x) = - t[(s 1 \qm i - 1 \P*)(s 2 \x(l - X -)\0) 

+ (s 1 \pm i \P*)(s 2 \x(l - x+)\0)} . (32) 

The action of rc_ and x + on |0) and of x on | — 1), |0), |1) is given by relations (17) and 
(18). The expectation of the ZRP-operator m k in the stationary state is calculated from 
the product measure as 

( Sl \m k \P*)=z (33) 

where the fugacity z can be identified with the probability for occupancy of a site in the 
stationary state. Equation (32) simplifies to 

(x) = - tqz(0 + 1) - tpz(0 - 1) 

= t{p - q)z (34) 

which in terms of the ZRP density c reads 

( x )=t(p-q) T ^-. (35) 

Taking the time derivative and substituting c by the £-ASEP density p, the average 
velocity of a tagged particle is obtained as: 

v=(p-9) i-(7-V (36) 
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3.4- Diffusive motion: Calculation of the diffusion constant D 

Equations (36) and (286) reduce the derivation of the diffusion constant D as defined 
in (21) to the calculation of the temporal derivatives of the first and second order term 
of the second moment of x. 

The first order term (x 2 )^ is obtained by a straightforward calculation similar to 
the one of (x)^: 

j t ( x Y ) = -2(p + q)z. (37) 

This is the direct contribution to the diffusion constant that one would have for 
Markovian dynamics of the tagged particle, i.e., in the absence of any memory effects 
resulting from the interaction. 

The derivation of the second order memory term of (x 2 ) as given in equation (27b) 
is a bit more subtle and requires an approximation to yield an explicit expression in 
terms of the stationary density p. Inserting (26), absorbing e - H ot^ e ~ H on j n ^ |p*^ anc [ 
into (s\ and taking the time derivative, it simplifies to 

^.(a;2)(2) = i. f'dn r dT 2 (s\x 2 Ve- H °^-^V\P*,0) 
dt dt Jo Jo 

= f* dT(s\x 2 Ve' HoT V\P*,0) . (38) 
Jo 

A straightforward calculation, where V is replaced by (24) and equation (12) and 
stationarity of the zero range distribution is used, turns the memory term (x 2 )^ into 
the integral of the current- current correlation function. Translational invariance of the 
system then yields the so far exact expression: 

4(^ 2 ) (2) = f &r 2 { Sl \ [gWi - 2pqm i +p 2 m i+1 ] e'^m^P*) . (39) 
dt Jo 1 J 

For symmetric hopping rates (v — -^{x) — 0) we define the diffusion constant by 

D = lim 4 (x 2 ) ■ (40) 
Applying the operator relation 

d t {si\ni = ^(si|[mi_i + m i+1 - 2m*] , (41) 

to expression (39), summing up (39) and (37), and without loss of generality setting 
v< 

(x 2 ) = (n (t)m (0)) P * - (n )p«(m )p* . (42) 



i = the time derivative of the second moment 4i{x 2 ) assumes the form 



dt 

d_ 
It 

This exact result is now approximated in terms of the stationary density c = (n ) by a 
linearization of the density at time t: 

no(t) =c + e (t). (43) 

The investigation of the time dependence of the fluctuations (e x (t)) in continuous space 
which is governed by the linearized operator equation (41) leads to a solution for 
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the time-dependent correlator in form of a Gaussian function for the initial condition 
(ej(O)mj(O)) = z. The evaluation at x = yields: 

(e o (0mo(0)> = z^p, = - JL= (44) 
V Ant v Ant 

In a last step, the stationary ZRP density c is replaced by the £-ASEP density p, and 
the time-dependent diffusion constant D is obtained as: 

D= LzlE 1 . (45) 

This is consistent with the well-known subdiffusive behaviour [14-16] as observed for the 
symmetric exclusion process (£ — 1). The novelty of our result is the density dependence 
of the amplitude. Treating the more general case of biased rates p ^ q one expects a 
diffusion constant which depends on the initial state chosen, see [20] for a recent similar 
result coming from a phenomenological approach. 

4. Hydrodynamic equation 

4-1. Derivation from ZRP properties 

The one-to-one mapping between ZRP and £-ASEP is now exploited in order to derive 
the hydrodynamic equation of the £-ASEP. In the ZRP there is no exclusion interaction 
between the particles and a hopping event occurs with hopping rates q or p, whenever a 
site is occupied by at least one particle. The probability to find site i non-empty at time 
t shall be called Zi(t). In the case of next-neighbour hopping and periodic boundary 
conditions, the ZRP density evolution at any site i is described by the master equation: 
d 

— Ci (t) = qzi-^t) + pz i+1 (t) - {p + q)Zi{t) . (46) 

In the hydrodynamic limit, the lattice constant a approaches zero on a coarse-grained 
scale. Substituting the discrete variable i by a continuous variable y, the master equation 
may be expanded into a Taylor series in powers of a: 

- aB —dy~ + a S ^y I ~ + ° {a ] (47) 

where 

B=p-q, S = P -p (48) 

characterize the contributions from the biased and symmetric part of the motion. 

Assuming local stationarity at sufficiently long times [17], the fugacity may be 
substituted by the stationary density, making use of the fugacity-density relation 
z(y,t) = i+c(y t } which is calculated from the well-established stationary properties of 
the ZRP. Inserting this expression into (47) yields a continuity equation for the ZRP 
density as a function of coordinate y and time t. The corresponding hydrodynamic 
equation for the £-ASEP is calculated by applying the transformation rules (16, 2b). 
The ZRP current J(c, d y c) which is a function of the ZRP density and density gradient 
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is thereby mapped onto the £-ASEP current j(p,d x p), being a function of the £-ASEP 
density and density gradient with respect to a time-dependent coordinate x(t). Carrying 
out all substitutions in (47) results in the sought-after hydrodynamic equation for the 
£-ASEP: 



dpjx, t) 
dt 



- aB 7T 

ox 



p(x,t)(i-e P (x,t)) 



+ a 2 S 



dx 2 



p(x,t) 



i-(e-i) P (x,t) 



l-(£-l)p(x,t) 

Equation (49) can be rewritten, introducing the 'effective density' 

P 



X 



l-(£-l)p 



and the hole density 

p h = \-ip 

as 



dp(x, t) d h 2 d 2 

~dT = - aB Yx [xp]+aS dx^ 



X 



(49) 



(50) 



(51) 



(52) 



d x j(x,p h )- 



The form (52) of (49) reminds of the hydrodynamic equation for the ASEP in the case 

e = v. 



(53) 



The concept of the effective density was first introduced in [18] where the average velocity 
of iV extended particles of length £ on L lattice sites is proposed to equal v — 1 — x- I 11 
terms of v, as derived in the previous section, equation (52) reads: 

d t p = - aBd x [vp] - a 2 Sd xx v (54) 
= - d x j(p,v) . 

Equations (52) and (54) match the result (36) for the average particle velocity v as 
obtained for a tagged particle in the previous section. 



4.2. Mapping between ASEP and £- ASEP 

A one-to-one mapping between £-ASEP and its special case of £ = 1 can be stated 
explicitly [4]. Let p 1 (x',t) denote the ASEP particle density as a function of time and 
coordinate x'. The corresponding transformations are 

x ' = x - (J - 1) [ X dxp{x, t) + 0(a) 
Jo 

p(x,t) 



p\x',t) 



X- 



(55 a) 
(556) 



l-(£-l)p(x,t) 

They prescribe a mapping between states on a lattice of L' sites containing iV particles 
of length 1 and a lattice of L = L' + (£ — 1)N sites containing the same number iV 
of particles which have length £. States which are transformed into each other will be 
called ZRP-equivalent states in the following. Scanning the lattice from left to right 
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and denoting each hole and each particle encountered in an ordered sequence 7$ with 
7i G {0,A}, where represents a hole and A represents a particle, yields identical 
sequences, whenever two states are ZRP-equivalent. 



5. S£/(2)-symmetry of the £-SEP 

For £ = 1 and in the case of symmetric hopping rates p — q — 1 the £-ASEP reduces 
to the symmetric exclusion process (SEP). One of the fundamental algebraic properties 
of the SEP is the fact that its Hamiltonian is S'L r (2)-symmetric [19]. In the following, 
the S'[/(2)-symmetry and its implications are investigated for the generalized symmetric 
exclusion process for particles of length £ (£-SEP). 



5.1. Quantum Hamiltonian formalism for the SEP 

In the case I = 1 the quantum Hamiltonian in terms of single-site particle creation and 
annihilation operators s k and number operators rik,Vk of the SEP is given by: 

,SEP T + ( 56 ) 

K,k+i = ~ s k s k+i ~ H-Sk+i + v k n k+ i + n k v k+1 . 
Choosing a tensor representation with the single-site basis 

n> = (J). |»> = (J) (57) 

for one-particle states \A) and empty states |0), representations of all operators s^ = 
(a x =f cr y )/2 and n k — (1 — o" 2 )/2, v — (1 + o~ z )/2 may be constructed of Pauli matrices 



a 



x,y,z 



The SU (2)-symmetry of the SEP has useful consequences: The overall particle creation 
and annihilation operators S T = J2k s t together with the operator S* 3 = X^(l/2 — n k ) 
form a Spin- 1/2 representation of the S'?7(2)-Algebra. 

Utilizing the fact that H SEP commutes in particular with S T , one can show that the 
local density pu = (n k ) = (s\nk\P(t)) satisfies a diffusion equation: 

j t (n h ) = -(n k H SEP ) = (n fc _!) + (n k+1 ) - 2(n k ) . (58) 

This relation reduces the density evolution of the many-particle problem to a single- 
particle problem (solution of the lattice diffusion equation) and also implies a 
correspondingly simple hydrodynamic limit, viz. the diffusion equation. More generally 
the SU(2)-symmetry implies that any fc-point many-particle correlation function can be 
calculated from an associated problem with at most k particles. 

Having asserted the existence of a general one-to-one mapping between ASEP and 
£-ASEP (556), naturally the question arises, which operator in the £-SEP state space 
follows satisfies the same kind of equation (58). 
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5.2. Quantum Hamiltonian formalism for the £-SEP 

The generalized form of the Hamiltonian (56) for i > 1 is: 
rre-SEP _ h e-SEP 

n — ' l k,k+i 

hk,k+£ P — - s k n k+\nk+2 ■ ■ ■ n k+ e-isl +e 

- sln k+l n k+2 ■ ■ ■ n k+i _is k+e (59) 
+ v k n k+ i ■ ■ ■ n k+e 
+ n k • • ■ n k+ i_iv k+ e . 

There are few formal but important physical differences in comparison with the SEP- 
Hamiltonian (56). For instance the action of s^ on any configuration of the lattice S is 
not equivalent to the creation and annihilation of an extended particle, but it creates 
or annihilates just a piece of it (one monomer). Furthermore H l ~ SEP is not anymore 
symmetric under the action of a true particle creation or annihilation operator and such 
an operator is of no relevance for the SU (2)-symmetry. A formula how to construct the 
operators which are taking over the role of S^' 3 for the £-SEP as a representation of the 
Lie-Algebra SU(2) is given in the following. 

5.3. Construction of a new 'creation' operator S~ 

The transformation between SEP and £-SEP relates states to each other which are ZRP- 
equivalent. 

Definition: Let \£' 2 ) be two arbitrary SEP states obeying 

\Q = s-\Q (60) 

and let |£ 2 ) denote the ZRP- equivalent i-SEP states. S~~ shall be defined by the 
relation 

I6> = 5-|&>- (61) 
S' + ' 3 are defined analogously. 

Due to their construction, S T ' 3 form a representation of SU(2) and the Hamiltonian 
is symmetric under their action: 

[H e - SEP ,S ± > 3 } = 0. (62) 

While S~ transforms between vectors of a fixed length 2 n , S~ maps a vector of length 
2™ to a vector of length 2 n+e ~ 1 . To be able to operate within one state space of fixed 
dimension, a new model of the SEP and £-SEP, generalized to two classes of particles 
A and B will be introduced. For the sake of simplicity of notation the case £ = 2 shall 
be considered first. The results are generalized to £ > 1 afterwards. 
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Figure 4. Creation of particles in a monomer and dimer system, yielding ZRP- 
equivalent states: S~ creates dimers on available and additional volume. 

The basis of the new state space K, is given by the tensor product states of the single-site 
basis 





(1) 




(o) 




(o\ 


|0) = 





, \A) = 


1 


, \B) = 















I 1 ; 



(63) 



The matrix representation of the single-site A- and 5-particle creation, annihilation and 
number operators in this basis is: 
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s B+ = 











s B - = 
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1° 
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1 














1 ) 



(64) 



Monomer space M. (£ — 1) 

The new representation of any monomer configuration G {^4,0} on the lat- 

tice 5" of length L' is realized by the following steps: 

• Double the lattice 5", adding L' new sites to its right end. 

• Fill the new sites with 5-particles. 

• Construct the tensor product state representation of the complete chain in the new 
basis (63). 

A monomer state which has previously been represented by a 2 L ' dimensional vector 
G H,® c \ is replaced by a 3 2i ' dimensional vector G /C. All new monomer states 
of arbitrary even dimension form a subspace of /C: 

= = G /C : a + v = b, £>-particles at right end} (65) 

where a, b and v denote the number of A-particles, -B-particles and holes respectively, 
which are contained in the configuration represented by 

The dynamics of the monomer system are governed by the monomer Hamiltonian H m . 
H m is obtained by substituting all H SEP operators in (56) by the ones which are labeled 
with a superindex A as introduced above. The new vacancy operator v is defined as 
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v = 1 — n A — n B . The action of the monomer Hamiltonian H m is local on the left half 
of the new monomer system and restricted to A-particles and vacancies. 



Dimer space V (£ = 2) 



The new representation of a dimer configuration e { A 0} on a lattice S of 

length L is constructed in the same way with the only difference that the number of 
lattice sites added equals the number v of zeros in {Pi}. Thus a state previously repre- 
sented by a 2 L dimensional vector |£) is replaced by a 3 L+V dimensional vector \5) G /C. 
The subspace of /C, consisting of all such dimer vectors \5) is determined by: 

^ = {($)} = {( K ) £ K. : v = b, A-particles exist pairwise, fi-particles at right end}. (66) 

The dimer Hamiltonian H d for a system of size L + v is given by the H 2 ~ SEP for a 
lattice of length L, where all operators except v are labeled with a superindex A. 



Mapping operators 



The new representation enables one to give the mapping between ZRP-equivalent 
monomer and dimer states explicitly in operator form: 

\5) = P%) (67a) 

|A*> = W (676) 



where 



p= n 

i=N-l 



l + (s 



B— „A+ 
N b N 



i+1 



Up, 



3-1 



l)n. 



j=N 



N-l 



p*= n 



i=i 



1 + ( II PjJ-i4~Sn + 
j=i+l 



(68) 



The permutation operator Pjj-i permutes two single-site vectors in the tensor product 
and can be expressed as: 



P lJ =v t v 3 +n A n A + n?nf + s A - S A + 



s A+ s A ~ 
b i b j 



Q B— Q B+ I „B+ „B— 
b i b j -T b i b j 



s A - 



,B+ B- „A+ 



^B- A+ A- B+ 
r>j bj bj 



(69) 



The construction of S = S A is now obvious: Let |<5j), (i = 1, 2) be ZRP-equivalent 
monomer and dimer states with 



Then: 



P t \^ 2 ) = P t S'\ f i 1 ) = P t S-p\S 1 ) 



(70) 
(71) 
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s 



Figure 5. Action of the mapping operators P, P* on dimer states (right) and monomer 
states (left), construction of the creation operator S~ from the monomer creation 
operator S~ in dimer space. 



5.4- Diffusion equation for i-SEP operators 

All dimer operators O 5 can be constructed as P t O tl P from monomer operators O^. The 
operator n k = P t n k P takes over the role of the monomer number operator: n k is the 
quantity which fulfills a diffusion equation with respect to the dimer Hamiltonian H d : 



The validity of PH d \5) = H m P\fi) in tC follows from the construction of the mapping. 
However, n k is not diagonal in the chosen basis of /C, and it is hard to draw conclusions 
as to its expectation value. Therefore a diagonal operator Q k is constructed, which 
equals n k in its action on (s\ = (s\P. 

The operator picks those dimer states from (s\, which are ZRP equivalent to the ones, 
rik picks from (s\. The action of an operator Qk which replaces hk thus cannot be local 
on the site k but must involve several lattice sites. Its action on a certain configuration 
depends on the number of particles and its label k. 

Qk can be expressed in terms of diagonal matrices only: 
fc-i 



{s\n k H d 



(s\PP t n k PH d 
(s\n k H m P 

(s|(n fc _i — njb+i + 2n k )P 
(s|P*(n fc _i - n k+l + 2n k )P 
<s|(nfc-i -nfc+i + 2n fc ) . 



(72) 



(73) 



a=0 
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where 



fc-i 



for a = : Q° k 



n v ° n k 



s=l 



for a > 1 : Q£ 



re^ Q \i=l / \sG5 r / 

• , r a ) : [n e {1, • • • , a + k - 2}] A [n + 2 < r i+1 }} 
k-l}\[Uf =1 { ri ,r l+1 }]. 




(74) 



n a = {r = ( n ,-- 

S r — {1, • • • , a + 



The local monomer number operator nk picks all states from (s\ where site k is occupied. 
Q k instead chooses all such dimer states from (s\ where the (a + l) th dimer, counted 
from the left, covers sites k + a and k + a + 1, and where a is a number between 1 and 
k. Thus each Q% sums up all possible configurations of placing a dimers (a A-particle- 
pairs) and k — 1 — a vacancies on the first k + a — 1 sites. The position of the left 
A-particle of dimer number % is chosen with the element of the vector r. Its elements 
must appear in certain configurations due to the pairwise arrangement and exclusion 
interaction of the A-particles. 

It is straightforward to generalize Qk to the case of particles of arbitrary length i > 1. 
The case of a > 1 in (74) is substituted by the general expression: 



6. Summary and Conclusions 

The purpose of this work was to investigate the properties of extended interacting par- 
ticles, moving stochastically on a one-dimensional lattice. The main results can be 
summarized as follows: One-to-one mappings between £-ASEP, ASEP and a certain 
class of ZRP have been stated explicitly. It has turned out very useful to exploit those 
transformations in order to derive basic properties of the £-ASEP, in particular the time- 
dependent diffusion constant for a tracer particle and the hydrodynamic equation for 
the local density evolution. A tagged £-ASEP particle shows the type of subdiffusive 
behaviour as known for the case £ = 1 of the ASEP. The extension of the particles as a 
new feature becomes manifest in the prefactor of the diffusion constant, which has been 
calculated as a function of the particle density. The methods applied are suitable for a 
generalization to any process which can be mapped onto some zero range process. As 
a most important outcome of the mapping, the hydrodynamic equation of the £-ASEP 
has been deduced from microscopic properties of the discrete system. The resulting 
nonlinear and convex current- density relation (for finite £) is qualitatively similar to 
that of the ASEP, but it also shows some new features: It is asymmetric for £ ^ 1. The 




(75) 



n a = {r = (r u ---,r a ): [n G {1, • • • , k + (£ - l)a - £}} A [n + i < r i+1 }} 
S r = {1, • • • , k + {I - l)a - £}\[Uf =1 {r u n + 1, • • • r i+{t -i)}] ■ 
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symmetry of particle and hole density is broken. The hydrodynamic equation has a 
natural form for £ > 1 if expressed in terms of the particle density and of a generalized 
average particle velocity. All results obtained for the £-ASEP also hold for a polydisperse 
system of particles of arbitrary length where the length parameter £ must be replaced 
by an average length £ = jq Yh=\ h- 

In the case of £ = 1, it is known how to link some hydrodynamic properties of the ASEP 
to the algebraic structure of the stochastic many-body system. Especially for the case of 
symmetric hopping rates (SEP), the S'L r (2)-symmetry has proved a valuable attribute. 
In this work, the S'[/(2)-symmetry has been established for the case of extended par- 
ticles (£-SEP). A formalism has been introduced in which all SEP-operators may be 
generalized to £-SEP operators under the condition of ZRP-equivalence. 

References 

[1] T.M. Liggett, Interacting Particle Systems (Springer, New York, 1985). 

[2] G.M. Schiitz, in Phase Transitions and Critical Phenomena 19, 1, 
C. Domb and J. Lebowitz (eds.), (Academic Press, London 2002). 

[3] T. Sasamoto and M. Wadati, J.Phys. A 31, 6057 (1998). 

[4] F.C. Alcaraz and R.Z. Bariev, Phys. Rev. E 60, 79 (1999). 

[5] G. Lakatos and T. Chou, J. Phys. A 36, 2027 (2003). 

[6] L.B. Shaw, R.K.P. Zia, K.H. Lcc, Phys. Rev. E 68, 021910 (2003). 

[7] L.B. Shaw, A.B. Kolomcisky, K.H. Lcc, J. Phys. A: Math. Gen. 37, 2105 (2004). 

[8] L.B. Shaw, J.P. Scthna, K.H. Lec, cond-mat/0403523 (2004). 

[9] C.T. MacDonald, J.H. Gibbs and A.C. Pipkin, Biopolymcrs 6, 1 (1968). 

[10] C.T. MacDonald and J.H. Gibbs, Biopolymcrs 7 707 (1969). 

[11] M.R. Evans, Braz. J. Phys. 30, 42 (2000). 

[12] F. Spitzcr, Adv. Math. 5, 246 (1970). 

[13] J. Buschle, P. Maass, W. Dietcrich, J. Stat. Phys. 99, 273 (2000). 

[14] S. Alexander and P. Pincus, Phys. Rev. B 18, 2011 (1978). 

[15] H. van Bcijcrcn, K.W. Kchr, R. Kutncr, Phys. Rev. B 28, 5711 (1983). 

[16] R. Arratia, Ann. Prob. 11, 362 (1983). 

[17] C. Kipnis and C. Landim, Scaling Limits of Interacting Particle Systems, 

(Springer, Berlin, 1999). 

[18] A. Ferreira and F. Alcaraz, Phys. Rev. E 65, 052102 (2002). 

[19] G.M. Schiitz and S. Sandow, Phys. Rev. E 49, 2726 (1994). 

[20] C. Godrcchc and J.M. Luck, J. Phys. A (2003). 



